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ABSTRACT 

The well known relationship between gauge invariance and current- con- 
servation is exhibited within the usual quantum mechanical formalism. It 
is then shown that the use of Gauge Invariant Atomic Orbitals does not 
necessarily lead to the expected current conservation. The reason is found 
to lie in the constrained nature of the gauge invariance which is provided 
by the use of GIAO's. It is concluded that this invariance is, of itself, 
no argument in favor of their use. 
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I. INTRODUCTION 

If 

If one views the Schroedinger equation from the point of view of 

i 1 

Lagrangian field theory, then there is a well known connection between 

gauge invariance and current conservation. In the next section of this 

paper we will establish a similar connection within the ordinary quantum 

mechanical formalism. In the third section we then show that the use of 

2 

gauge invariant atomic orbitals (GIAO's) does not necessarily lead to 
the expected current conservation, and the reason is discussed. This 
then leads us in the fourth section to examine in detail the nature of 
the gauge invariance which is provided by the use of GIAO’s. We conclude 
that, contrary to what is probably the generally accepted view, this in- 
variance is, of itself, no argument in favor of their use. The body of 
the paper is concerned with static fields and with energy eigenfunctions. 

An appendix is then devoted to a discussion of time dependent situations. 

II. GAUGE INVARIANCE AND CURRENT CONSERVATION 

We consider the problem of N non- relativistic electrons moving in 
external static electric and magnetic fields, and we confine attention 
to energy eigenstates. We now make use of the following general criterion 
of invariance: Since the formalism of quantum mechanics is invariant to 

unitary transformations of wave functions, one will have invariance to 
particular transformations if these can be shown to lead to unitary trans- 
formations. It is well known that gauge transformations meet this criterion. 
Namely if the vector potential undergoes a gauge transformation 



where is some real function, then the wave function undergoes the 

unitary transformation 


A A A 

C a') = e * c t) 


u) 


and correspondingly the Hamiltonian also undergoes a unitary transfor- 
mation 

* .4 <0. 


v+ (a 1 ) ^ a VtCt) e 


( 2 ) 


where (we use atomic units and c is the velocity of light) 


A. » 12 U (3) 

Co 

In particular we will be concerned with families of gauge transformations 


X C V ? •=. 
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where the ^ are arbitrary real numbers, and where the "X o( are 
some given real functions. 

If now we expand V4C/^7 ih powers of the ^oi we find through 
first order in the 


l-\ ( a‘> C \ Uffol * 


(5) 


where [ , ] denotes the commutator. Then from first order perturbation 
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where 
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and where 

A* - t %- ^ 

c <-*1 

However since (2) is a unitary transformation of H - , the energy is 

C A) 

gauge invariant, whence must vanish. Therefore we must have 

t A t, Ao4^ 14 ^ *=-0 (?) 

which is of course true since \4f C A*.} ^ ^ • Thus 

3 

the hypervirial theorem (7) can be viewed as being a consequence of 
gauge invariance. 

It is now easy to show, after some integrations by parts, that (7) 
can be written as 

'"■e 

S 3 * =0 ( 8 ) 


where 

all 


j is the one-electron current. 

X c* we have 


"= 5.0 


Then since (8) is to hold for 


(9) 


i.e., current conservation. 

So far all our discussion has concerned exact wave functions. We 

now quote two relevant results for variational wave functions: 

(i) If the set of trial functions is invariant to the set of gauge 
-» v- J\- 5 

transformations & then the optimal trial functions will 

satisfy the hypervirial theorems (7) and hence (8). 
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(ii) We will have gauge invariance in the sense that if we use this 

same set of trial functions with both and C4i/^ , 0 . then^ we 

y ms 

will get the same energies and also the optimal trial functions will be 
related by formulae like (1). 

si 'v'S/cr- c 7 ■■■ ic> c.tr " ■ -tj .? & a \ >. } s%n.:A uavawo.h 

In particular it is easy to show^ that unrestricted. coupled Hartree- 

j.:-; 28i.r,'v ‘.--.r .-lair.-'v, jtif ..j aoaadw 

Fock theory meets these conditions for all , and indeed in this 

/•f _ i - ; tc T f 5 .’ : ‘ .{ . f *. 

theory the energies - are gauge invariant ,- ; arid ;, the t one -electron current i£ 
8 

conserved.: .However in mos.t, variational calculations this will .not be 

f , \ ' ■ e t>-.X 3 *mjy 3 in 


the case. 

; J: - >3 a ;■?. p c j ar i* a \ » & • v i v J 
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III. GIAO'S AND CURRENT CONSERVATION 




t '■ O 


# * o.* vc-d'ij v.*i s-ih..=c: r-;fs> «>:• 

For clarity we will consider a very simple system: a single electron 

tf.jc . 1 1 ’V U j' 1 ,r -if’ 

moving in the field of two fixed centers, and in a uniform magnetic 

<V"-' V ^ 

field B , Further, in accord with what seems to be the custom in quantum 

chemistry, we will confine attention to a very special class:; of gauge 

cl^i M si •'•-■) oral*. an • ai •. 

transformations. Namely for the vector potential we take / 

v. ‘ - 


A = 




(10) 


and then confine attention to the gauge transformations 

ft oi.' •-' ! v . c;. - f 

%! . ■ 5 • 

> eva-.- -rr> 'm. : ' s - *. . > 

where d is an arbitrary constant vector. From (11) we have 

$3 \ a; * Y c-n.‘ t ? -• x 1 — i... ^ - 

, . . . | , V 

:.v f noilo/’ ■ 'v^-i i.M:* qiv ^ - * 

A i ” 


Ui) 


( 12 ) 
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Now since the use of GIAO’s is known to give gauge invariant energies 
(we discuss this further in the next section) one might suppose, from 
our earlier discussions, that one would necessarily have 


(VS v* -o 


(17) 


..-•r and :• • 


" h f A J so ’ ' 

We will now show that this is not the case. 

First we consider (18) . One readily finds that 


) ^ A ) - ^ ) k ) 4 M ' X f 4^ (frj- |B ) ^ -j 


(18) 


-r 




^BR 


f / d (<$>,. e 2 




/pg+^] . ^) 4 c, c7\ (19) 


Our assumption that 4? and 


are s-states then makes the first 


line of (19) vanish, and allows the second line to be written as 
We now note that for B = 0 the expression in the brackets becomes 


— ( Vk 

2 ^ 


) 


which in general doesn’t vanish. Therefore . (19) will in general vanish 
only if 0 ( is real. 
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We now turn attention to ^ • Then assuming 

real and invoking the symmetry of and one readily finds 

that 



and in order for this to vanish the bracket must vanish (unless $ 
which seems unlikely. In particular if we expand in powers of B the 
leading term in the bracket is 


, ZC . 




( 20 ) 


■3>v. • 


.and there seems no reason for even this to vanish, and indeed $ calculation 
with Gauss ians 

. _ A h y v >»*-0 

yields a non vanishing, result (unless ^ ) . Therefore we conclude 

that in this particular example we cannot, have both (17) and (18) , and 
hence we conclude that in general the use of GIAO's does not guarantee 
such theorems. 

The reason for this failure is not hard to find. As stated in (i) 
of Section. II, the theorems will hold if the set of trial functions is 
invariant to the appropriate gauge transformation. However our' 1 set, 
which (for fixed ^ and ) consists of the single function 

, evidently changes under the gauge transformation (12) : Namely 


from (15) we have 



- / i. * Vt* A v 

^ ^ ^ ^ C. (> ^ sp 




( 21 ) 


Thus the conditions of (i) are not met and so phe theorems don't hold 
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Another, equivalent, reason for failure, one which makes direct 
contact with (6) (the argument pb^ve of ppurse dealt directly with (7)), 
is the following: First order perturbation theory in the fora pf equa- 

tion (6), is a special case of the generalized Hellmann-rFeynman theorem 
for the parameter . Now one knows*^ that this theorem vfll also 

hold in variational calculations provided that the set of trial functions 
is independent of the . Evidently this is not the case here. 

Thus the first order energy change is not given by a formula like (6) 
but rather by one which contains extra terms, and therefore although fhe 
first ordpr energy change does vanish fn accord with gauge invariance* 
that vanishing dees net lead to (7), 

These considerations naturally raise the question* what then is the 
nature of the gauge invariance which is provided by the use of (JIAO' s? 

It is thfe question that we now turn. 

IV. GIAO'S AND GAUGE INVARIANCE 

.ft 

Clearly the fact that t changes with A , the fact which led 

to the failure in the previous section, is precisely the feature of 

A. 

which leads to the gauge invariance of the energy. Thus the 
transformation (?1) is exactly the transformation (1) (for a one electron 
problem) and hence one hap 






i.e. , the energy jls gauge invariant. However having recognized that the 

gauge invariance is guaranteed by the way in which the orbitals change 
14 

with gauge, we then note that by a similar device we can make any 
calculation equally gauge invariant. Thus for example consider instead 
of the function 


t 


h 


-* tlfo +■ 


-i teg:** 


- e 


where 


t 


and 




I emu i ^ are any given vectors. Then since we have 

_ . i -£■ X r a 
^ t a V 


a calculation based on ^ will also yield a gauge invariant energy. 

In particular if we choose = 0 then . evidently we could con- 

clude that even the use of simply 

♦‘I*'* W - 

with A = •^'53*.* and e ° C^, {v) otherwise, will yield gauge 

invariant results. 

More generally it is clearly always possible to get such, what we 
might call "enforced gauge invariance" or "constrained gauge invariance", 
for the energy calculated using any set of trial functions whatsoever in 
any one gauge simply by changing the trial function appropriately as one 
changes gauges (or equivalently, by always returning to the original 
gauge to do the calculation) . Further it would seem that if one proceeds 
in this way, and as we have seen in using GIAO's one does precisely this,^ 
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that then the resulting gauge invariance should not be thought to lend 
any special credance to tfye calculation since requiring it has simply led 
to a specified connection between the sets of trial functions to be used 
in different gauges. Rather one's a priori, expectations concerning the 
accuracy of the calculation, i. e. , one s a priori expectations concerning 
tlie physical completeness of the set in any one gauge, must clearly be 
based on other considerations . In contrast, if in the spirit of (i) and 
(ii) of Section II, the set; of : trial functions is invariant to gauge 


transformation^ then equally clearly this is a relevant piece of informa- 
tion since evidently it does supply some information concerning the 
completeness of the set. 

Now in fact, as is well known, there is a^real physical basis for 

is, through 




plus uniform magnetic field problem which for B = 0 has 
16 


the use of GIAO's: The function 

first order in B , the exact eigenfunction of a one-electron one-center 

^ L as an 

eigenfunction. Moreover we would emphasize that nothing we have said 
about "constrained" or "enforced" gauge invariance should obscure th>5 
point. Rather we are saying that this, and not the resulting gauge 
invariance, is what suggests that the use of GIAO's may be an effective 
tool in dealing with problems of molecules in uniform magnetic fields. 

On the other hand, of course, the failures which we have discussed are 
also real and should be weighed in the judgement. 


Returning to our special, example, since any choice of 


and 


will yield a gauge invariant energy 'which however depends on 
and * it .is natural to ask for the best ^ and 


y%* 


in the 



XX 


sense of the variation principie. To find out we of course simply let 


$ 

% 


iTP . 

and be variationai parameters. If is then of interest 


to note that the resultant set of trial functions is now invariant to 


gauge transformations since y 


A' 


can be gotten from "V 


X>y 


replacing J® and 5*- by and d respectively, 

so that a gauge transformation simply turns one trial function into 
another. Thus from (i) and (ii) of Section II it follows that by intro-? 

ducing this extra freedom"^ we will ensure that (17) an$ (18) are now 

'Of-' 

and fV 

are not (0,0,0) and (0,0, R)) , and also ye will still have gauge} invariant 
energies and properly gauge variant wave functions, but now in Sty un- 
forced way (one uses the same set of trial functions in all gauge?) rather 
than by constraint. 


V 1 f ' . i 

satisfied (evidently therefore in our example the optimal 
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APPENDIX 


Although in the body of the paper we are concerned only with station^ 
ary states it would seem useful to add some comments on the time dependent 
case. Here under a gauge transformation the vector and scalar potentials 


change according to 
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!f (A_1> 

x c "3t 

Again invariance is expressed by the fact that the wave function undergoes 
a unitary transformation 

-W£\ $ ( ) ^ t (if, §) < A - 2 > 

and correspondingly 

■ ; H-C &-*’*’ Xtlty 5)e vA +i^A : (A-3) 

c i)V 

where W ^ 

A ^ 2- >C 

u-V 

If now we choose 

then the analogue of the result (6) for the first order energy change is 
provided by the time dependent Hellmann-Feynman theorem which yields 
a relation between the change in and the change in H , due to a 

change in a parameter, whatever the nature of that change may be. As we 
have emphasized, gauge invariance is expressed by the particular nature 
of . the changes foo) and in particular by the fact that(j&Ht^ is a unitary 

transformation. For our particular case then the general theorem 
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where q-' is a real parameter yields 

L% IV*) = ( 4 -, f) CH ; Ajf) (a-4) 

This result, a special case of the generalized Ehrenfest's theorem or 
time dependent hypervirial theorem, is of course a well known consequence 
of the Schroedinger equation. However we now see that it can also be 
regarded as a consequence of gauge invariance, and evidently in the time 
w\ dependent case, and with the £<* = 0 , it reduces to (7). Further it 
is easy to show that (A-4) can be rewritten as 


r&jjK =o < A - 5 > 

where ^ is the one electron density. Then since (A-5) is true for 
all 'Xet we have 

-o 

at 

i.e., current conservation. 

All this is for exact wave functions. For variational wave functions 
19 

one can show in complete analogy to the static case, that if the set 
of trial functions is invariant to the set of gauge transformations then 
the optimal trial functions will satisfy (A-4) . This condition is met 
by unrestricted coupled Hartree-Fock theory for all ^ and hence 


one has current conservation in that theory. 
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